Patterson Function from Low-Energy Electron Diffraction Measured Intensities and 

Structural Discrimination 
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Surface Patterson Functions have been derived by direct inversion of experimental Low-Energy 
Electron Diffraction I-V spectra measured at multiple incident angles. The direct inversion is com- 
putationally simple and can be used to discriminate between different structural models, (lxl) 
YSi2 epitaxial layers grown on Si(lll) have been used to illustrate the analysis. We introduce a 
suitable R-factor for the Patterson Function to make the structural discrimination as objective as 
possible. From six competing models needed to complete the geometrical search, four could easily 
be discarded, achieving a very significant and useful reduction in the parameter space to be explored 
by standard dynamical LEED methods. The amount and quality of data needed for this analysis is 
discussed. 
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Low-Energy Electron Diffraction (LEED) takes advan- 
tage of strong interactions between LEED electrons and 
atoms to achieve highly sensitive and accurate surface 
structure determinations [[!]. The prize to pay for that 
sensitivity based on a strong interaction is that simple 
first-order perturbation techniques are not enterily ap- 
propriated, and multiple-scattering is unavoidable to an- 
alyze the experimental data. Because multiple-scattering 
implies self-consistency Q, direct inversion of data to ob- 
tain the structure is quite difficult, and the most popular 
route to surface structural determination relies on a trial- 
and-error procedure, that becomes too demanding when 
the parameter space to be explored is too big, or the 
structure cannot easily be guessed beforehand. There- 
fore, there is a genuine interest in developping and explor- 
ing direct inversion procedures, regardless of small prac- 
tical difficulties or some limitations that might be found 
associated with them||. For LEED, an already succesful 
inversion technique based in ideas derived from holog- 
raphy exists || ||. An useful alternative to holographic 
ideas, based in the Patterson Function (PF) concept, has 
also been proposed in the literature by Tong et al. ||. 
These two proposals are based in different physical con- 
cepts, but have in common that they yield useful infor- 
mation on the surface structure by directly performing 
on the experimental data a Fourier transform-like with 
an appropriate kernel. The main desired effect of this 
mathematical operation is to filter out high order scatter- 
ing terms. Out of necessity, this is only an approximate 
procedure, but previous experience with holography, it- 
erative methods based in maximum entropy 0, or even 
quasi-dynamical methods || have shown that multiple- 
scattering is not necessarily a roadblock. 

In this paper we explore the PF idea, and by consider- 
ing a practical example based in experimental LEED-IV 
data measured at different energies and incident angles 
we show how useful these ideas can become when they 
are used to complete a real surface structure determi- 



nation. The strategy to defeat the limitations imposed 
by multiple-scattering is to collect data simultenously 
at different energies and incident angles. This results 
in overdetermination of the inverted data, where coin- 
cidences in features not directly related to the under- 
laying geometrical structure decrease as the data base 
increase Therefore, the effort to get the structure 
shifts from the trial-and-error search to the measurement 
part, which is not a real burden taking into account mod- 
ern LEED data acquisition techniques. 

The surface Patterson function (PF) can be com- 
puted by an appropriate phase sum of the measured 
intensities Hi: 



ki 3\\ 



(1) 



where ki is the wavevector corresponding to incoming 
electrons, q represents the momentum transfer for outgo- 
ing electrons detected in the LEED screen [kf = ki + q), 
and gu labels the different beams. Simple real-space scat- 
tering models analyzed by Tong et al.[^| help us to un- 
derstand why the multi-energy and multi-angle approach 
can reduce spurious features appearing in the stationary 
phase condition implied in Eq. (|l|). The PF yields max- 
ima at positions related to interatomic distances: we shall 
find that these can be used to discriminate reasonable 
from unreasonable structural candidates. 

To better understand how useful these ideas might be 
for real structural determination we take advantage of 
the experimental LEED I(V) curves we have recently 
measured at normal and non-normal incidence for (lxl) 
YS12 layers grown on Si(lll) substrates. These mea- 
surements have been used to perform a standard I(V) 
study of the structure, where the best correlation was 
found at a Pendry R-factor pf| of Rp = 0.21, which is 
accepted to be a low-enough value to credit the struc- 
tural model. Details related to this model have been 
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published elsewhere |Q, but the structural model is pre- 
sented in Fig. (|l|A) for reference. Basically, the surface is 
made of Si(lll) planes 180° rotated around the Si(lll) 
surface normal, with Y atoms located below the surface 
on T4 sites (w.r.t the Si(lll) substrate). The small re- 
laxations between layers included in the model are not 
crucial for the analysis presented in this paper, because 
the spots on the PF already have an intrinsic width on 
the same order or bigger than the values associated to 
the pure geometrical surface relaxations. 

Experimental LEED I(V) curves were recorded at nor- 
mal and non-normal incidence. The latter was achieved 
by moving the surface normal in a plane defined by the 
surface normal and the [121] surface direction. Mea- 
surements at 9 = 0°,4°,9° and 14° where taken for the 
same azimuthal angle, <f>. 19 beams (6 of which are not 
symmetry-related) where measured from 50 to 450 eV, 
except for normal incidence where only 18 beams could 
be collected (the specular beam was blocked by the elec- 
tron gun). Therefore, the total energy range for inde- 
pendent data is more than 8000 eV. The experiment was 
performed at room temperature: in the standard LEED 
analysis it was found that fine tuning of isotropic vibra- 
tions associated to different layers could make a strong 
impact on the Rp value, decreasing its value from around 
0.4 for r.m.s. displacements taken from bulk values, to 
around 0.2 for the final optimized surface valuesjllj. At 
the moment these effects are not taken into account in 
Eq. ([[]) , but we anticipate their main effect is to broaden 
the PF increasing the values of the correlation factor we 
have found below. 

The PF gives the interatomic distances between all 
atomic pairs in the structure. To simplify the discus- 
sion and make a better representation of our results, we 
have chosen the [10 1] cross-sectional plane, but the same 
conclusions would be obtained using a different plane. 
Fig. (QB) shows the interatomic distances for one unit 
cell calculated from the structure in Fig. (EK): stars cor- 
respond to positions where a perfect stationary phase 
condition would produce spots in the PF from I-V spec- 
tra. The corresponding atomic pairs are indicated next 
to the stars, using the same labelling as in Fig. (0A). 
This is a mere geometrical PF, and it is the best result 
that could ever be recovered for a given structure. There- 
fore, we compute correlations between experimental PF 
obtained from Eq. (Q) using all measured diffracted inten- 
sities and geometrical PF corresponding to different trial 
structures. As shown in Fig. (||), the experimental PF is 
made of spots with a certain width derived from factors 
like the limited domain of integration, atomic vibrations, 
defects, etc. 

A visual comparison between the stars in Fig. (^B) and 
the spots in Fig. (jfy induces us to introduce a correlation 
factor between PF's to facilitate their comparison. The 
new correlation factor might be defined in different ways, 
but to demonstrate the concept we shall here introduce 
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FIG. 1: (A) Upper pannel: Schematic positions and labelling 
of atoms in YSi2/Si(lll). The geometry corresponds to the 
best fit determined to experimental LEED data measured on 
the same systemjnj (B) Lower pannel gives the geometrical 
Patterson Function related to structure in upper pannel. 
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FIG. 2: PF from Eq. ([lj, and the experimental I(V) curves 
measured for 19 beams at four different incident angles. 



the simplest one, to keep the discussion as independent 
as possible from particular details of this definition: 



Rpf — 



S. 



spots 



ED 



Spf N, 



(2) 
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where D is the distance between stars corresponding to 
a given structure and the center of the closest experi- 
mental spot around it, N geo is the number of interatomic 
distances considered, S spo t s is the area covered by the 
spots, and Spf is the total area considered for the anal- 



ysis. The factor 



gives the probability to find a spot 



ZD 

N aco 



measures 



in the considered area, whilst the factor 
how far are the experimental spots from the correspond- 
ing stars. In order to preserve the meaning of this partic- 



ular Rpf, the factor 



Spi 



should not be too close to 1 



(e.g., spots covering the whole surface) or to (e.g., too 
few spots). By constructing model geometries at ran- 
dom, we estimate typical values for Rpf to be around 
0.25-0.20 when the structures are not correlated. 

The structural analysis for YSi2 epitaxially grown on 
Si(lll) requires considering six different models for dif- 
ferent stackings of the atomic layers. In our previous 
analysis of this structure) 11 , a detailed Rp analysis cov- 
ering different values for all the considered parameters 
has been performed. This extensive search reflects the 
sensitivity and highly non-linear dependence of inten- 
sities with structural parameters in multiple-scattering 
calculations Jl2[ . This is not the case for PF, where a 
stationary phase condition based in two first-order scat- 
tering events is sought. Therefore, the advantage of com- 
puting correlations between experimental PF and differ- 
ent models is double: (i) there is no need to run costly 
multiple-scattering calculations for the whole parameter 
space, and (ii) because the analysis is based in first- 



order scattering events, the Rpf is expected to behave 
smoother as a function of the involved parameters, and 
less local minima are expected. Although, on the other 
hand, a smoother variation may result in less sensitivity 
to different structures. Correlations between the exper- 
imental PF and the different models considered in the 
previous structural analysis are given in table |. To com- 
plement Rpf, we give next to it a number counting the 
percentage of coincidences between the stars and the PF 
spots. This has been defined as the number of stars which 
are inside the region in black in the figures. This number 
can be helpful to discriminate between models, as can 
been seen in Table |, but it has not been included in the 
Rpf definition because it depends on the width of the 
spots (a function of several factors, e.g. the pure thermal 
vibration of atoms is expected to contribute an intrinsic 
width in the order of several tenths of Angstroms). 

From entry named Experiment in table || we can in- 
mediately reduce the available six models to two (i.e., a 
66% saving has already been achieved by simply Fourier- 
transforming the experimental data). Differences in Rpf 
between model B and F, however, are not big enough to 
discriminate between them. In fact, the lowest Rpf cor- 
responds to model F, whilst the best fit based in Rp 
corresponds to the second lowest i?pp-value, model B. 
Notice that in both cases, the percentage of coincidences 
are among the highest. 

To understand the origin of this different ordering in 
models based in PF analysis, we test the consistency in 
the definition of Rpf by replacing the experimental I(V) 
curves by theoretical ones corresponding to the structure 
B. Results shown in table | under the column named 
Theory (19 beams) show that the definition for the R- 
factor is internally consistent and now the best model is 
clearly the correct one. Therefore, the question arises of 
what is the effect on the PF of exchanging sets of I(V) 
curves (i.e., experiment and theory) that differ them- 
selves by Rp = 0.21. We notice that multiple-scattering 
cannot be the source for the discrepancy as both, the ex- 
perimental and the theoretical intensities include all the 
multiple-scattering relevant for the problem: more likely 
this is due to residual factors (like anisotropic vibrations 
of atoms near the surface, defects, etc) or the size of the 
integration domain itself. Therefore, we turn to a visual 
comparison of Fig. (0) and Fig. (0A). Both PF's have a 



similar probability to find a spot: = 0.40 for exper- 



iment and 



0.39 for theory, but widths of individ- 



ual spots are too wide to obtain an accurate correlation 
with stars related to models B and F. Due to spurious 
features, both models become undistinguishable: this is 
not too surprising as the two models in question are very 
similar. 

Two simple ways to narrow these widths come to mind. 
From an experimental point of view, measurements made 
at low temperatures or recorded on samples as perfect as 
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TABLE I: Rpf and coincidences for different stacking models. Results are shown for the experimental data and two theoretical 
simulations (19 and 91 beams). Coincidences give the percentage of geometrical stars inside experimental or theoretical PF- 
spots. 



Code 


Structural Model 




Experiment 


Theory 


[19 beams) 


Theory (91 beams) 






Rpf 


(% Coinc.) 


Rpf 


(% Coinc.) 


Rpf 


(% Coinc.) 


A 


T 4 Rotated 0° 


0.159 


40 


0.114 


60 


0.113 


67 


B 


T 4 Rotated 180° 


0.140 


50 


0.069 


88 


0.026 


88 


C 


H 3 Rotated 0° 


0.175 


40 


0.122 


60 


0.146 


53 


D 


H 3 Rotated 0° 


0.150 


47 


0.123 


59 


0.060 


65 


E 


Atop Rotated 0° 


0.161 


50 


0.088 


73 


0.107 


86 


F 


Atop Rotated 180° 


0.136 


67 


0.083 


73 


0.030 


80 


G 


random model 


0.250 


13 


0.250 


13 


0.180 


6 




[121] direction (A) 



FIG. 3: (A) Upper pannel shows the PF obtained from 
Eq. (|l|), and the theoretical I(V) curves for 19 beams at four 
different incident angles computed for geometry in Fig. (|l|). 
(B) Lower pannel shows the PF obtained from Eq. (hi), and 
the theoretical I(V) curves measured for 91 beams at four 
different incident angles. 



possible are needed. On the other hand, the station- 
ary phase condition implied in Eq. [l] requires covering 
as much volume as possible in the Fourier-transform do- 
main (i.e., measuring as many beams at different angles 
for as many energies as possible). 

Our experimental setup puts some limits on the num- 
ber of beams we can measure, and restrict us to room 
temperature. We comment in passing that analyzing a 
lxl put us in the worst possible case for PF analy- 
sis. For the same energy range, any other superstructure 
will provide extra beams to be measured. Furthermore, 
the lxl usually would contribute less stars than a sys- 
tem with a surface superstructure, where different pair 
distances are involved. Finally, beams in the lxl mix 
surface and bulk information, a disadvantage that the 
fractional beams corresponding to a genuine surface su- 
perstructure do not present. 

Therefore, we explore the second requirement by sub- 
stituting our experimental I(V) curves by theoretical I(V) 
curves computed for model B for all the possible exit- 
ing beams (i.e., 91 beams in total). Indeed, comparison 
between Figs. (|| A) and ^ B) demostrates that the 
widths of the PF spots narrow when the data base is in- 
creased: the probability to find a spot has been reduced 
to s " pots = 0.25, nearly a factor of two w.r.t. the the- 
ory or the experiment based in 19 beams. Clearly, if we 
can measure more beams we could obtain more accurate 
correlation between stars and PF-spots. 

In conclusion, we have found that the PF concept al- 
lows us to significatively cut the parameter space to be 
explored in a real structure determination based in LEED 
data measured at different angles and energies. The extra 
burden of having to measure more beams and the require- 
ment for liquid Nitrogen cooling temperature should not 
be excesive for modern LEED equipment. The advantage 
is significative. The analysis based in the PF is very fast 
compared with the standard procedure, as it mainly im- 
plies Fourier-transforming the measured data. Inside the 
restricted parameter space determined by the PF analy- 
sis, it seems appropriate to run multiple-scattering cal- 
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culations. This is certainly more costly than the simple 
integral involved in PF but it might be worthwhile: its 
intrinsic non-linear dependence with the structure, and 
the presence of higher order correlations, will help to ob- 
tain accurate and trusty values for the final parameters 
of the model. 

Work supported by the Spanish Ministry of Science 
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